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We determine the total state dynamics of a dephasing open quantum system using the standard 
environment of harmonic oscillators. Of particular interest are random unitary approaches to the 
same reduced dynamics and system-environment correlations in the full model. Concentrating 
on a model with an at times negative dephasing rate, the issue of "non-Markovianity" will also 
be addressed with the emphasis on information obtained from the dynamics of the total state of 
system and environment: making use of criteria that allow us to distinguish between classically 
correlated and entangled total states, we employ a simple measure for the correlations emerging 
from the increase of the two local entropies, and relate it the nature of the correlations. 
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I. INTRODUCTION 

In a realistic description of open quantum system dy- 
namics one has to account for influences of the external 
environment [iHll • Despite the unitary time-evolution of 
the total state of system plus environment, the dynamics 
of the system itself will in general be non-unitary. Grow- 
ing correlations between the system of interest and its 
surroundings lead to a decay of the initially present co- 
herences. This line of thought is at the heart of decoher- 
ence theory and is put forward to explain the appearance 
of classical properties in quantum systems Deco- 
herence in particular is of relevance for quantum tech- 
nologies trying to make use of the vast computational 
potential forecast to applied quantum information pro- 
cessing 0. 

In the regime of weak system-environment coupling 
and short environmental correlation times the dynamics 
of an open quantum system may be described in terms 
of the Born-Markov approximation. The corresponding 
Markov master equation then has a generator in Lind- 
blad form [l^, [HI . The time-evolution of the system de- 
pends on its present state only. Often, however, such an 
approximation is not justified. Then, memory effects — 
typically incorporated by means of integrals over the past 
[S [13 — start to play an essential role. Yet, it is known 
that for arbitrary bath correlation functions exact time- 
local master equations can be derived 0. An impor- 
tant example is given by the exact master equation for a 
damped harmonic oscillator bilinearly coupled to a bath 
of harmonic oscillators [13, [HI • Here, we focus on a de- 
phasing open quantum system using the standard har- 
monic oscillator environment. 

In more recent developments, the question of how 
to distinguish Markovian from non-Markovian dynam- 
ics from a local (system) perspective was addressed. The 
analysis has been based both on a single snapshot of the 
dynamics and on the full time evolution of the open 
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quantum system within a certain time interval [13, [3 ■ 
In the latter approach, memory effects associated with 
non-Markovian dynamics are expected to cause tempo- 
rary increase in the distinguishability of states (in terms 
of trace distance, e.g. Under Markovian dynamics, 

on the other hand, the decay of distinguishability will be 
monotonic throughout. 

The role of system-environment correlations in open 
quantum system dynamics has raised some notable in- 
terest lately. In the context of quantum discord fl9| . 
for example, total states with no quantum correlation 
(zero discord) were shown to be the most general class 
of initial states allowing for completely positive reduced 
dynamics [l^, [2l| . Another interesting point is the re- 
lation between open quantum system decoherence and 
system-environment entanglement. Here, it is possible 
that the system decoheres completely without becoming 
entangled with its environment at all [l^, [2^ . 

Talking about the dynamics of open quantum systems, 
the flow of quantum information from system to envi- 
ronment and vice versa is an often employed, and cer- 
tainly intuitive picture. However, since there is a lack 
of a quantitative description, e.g. in terms of a continu- 
ity equation, such a picture should be used with caution. 
Care has also to be taken with respect to the need for a 
proper environment. At this point it should be noted that 
for single-qubit decoherence the dynamics may always be 
described in terms of stochastic fluctuations of external 
fields, i.e., the dynamics has a random unitary repre- 
sentation [U, [2^ . The dynamics may thus be modeled 
without invoking a quantum environment at all. Higher 
dimensions are needed to see proper quantum dynamics 

In order to shed light on the role of "quantuni- 
ness" in open system dynamics in a non-trivial (non- 
Lindblad) regime, we here investigate single-qubit de- 
phasing |23l.l27ll28|. by means of two different aproaches: 
we use a proper quantum model and also show how to 
find a random unitary representation of the same reduced 
dynamics. In the quantum case, we employ a useful rep- 
resentation of the total state of system plus environment 
discussed in [2^. This enables us to study the dynamics 
of both system and environment purity, as well as the 
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dynamics of system-environment correlations, which we 
also relate to non-Markovianity. 



II. QUANTUM DECOHERENCE MODEL: 
REDUCED AND FULL DYNAMICS 



Continuing our earlier work on the dynamics of system- 
environment correlations for a dephasing qubit |23| . we 
start from a typical model [2^,|3^ with total Hamiltonian 
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by coupling a qubit non-dissipatively to a bath of har- 
monic oscillators through the choices (sij - fssj 
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Here Q denotes the energy difference between the qubit 
states and the coefficients gx describe the coupling 
strengths between the qubit and each environmental 
mode of frequency ujx and annihilation and creation op- 
erators ax,a\. 

As environmental initial state we assume for 
each oscillator a thermal state p^horm = ('^^ + 
1) exp[—hu!xa\ax/kBT] with the mean thermal occupa- 
tion number nx = {exp{hu!x/ kBT) — at temperature 
T. Initially, we assume no system-environment correla- 
tions such that the total initial state is simply given by 
the product ptot(O) = Psys <8) Ptherm- Accordingly, for the 
reduced system state Prcd(i) = Tronv [ptot (01 the dynam- 
ical map £{t,0): Prod(O) Prod(i) is completely positive 

with Prcd(O) = Psys- 

Already at this stage we emphasise that on the re- 
duced level, this dynamics can equally well be described 
by random unitary dynamics as will be elaborated upon 

in Sec. nni 

The above quantum dephasing model ([Ij may be 
solved without any approximation. A possible approach 
to the time-local master equation is provided by the 
non-Markovian state diffusion approach to open sys- 
tems [1^ [S^l ■ We find for the reduced density operator 

Prcd(0 = Tre„v[ptot(i)] 

Prcd = Y Fz, ProdJ Y' 

This equation is solved by 
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V{t) = exp 
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and where the py represent the initial state of the qubit. 

Equations ([3]) and ^ involve the time dependent de- 
phasing rate "f{t) which by means of the spectral density 
of the environment J(w) — J2x=q I.9aP'5(w — lox) can be 
written as 



lit) 



ds 



dujJ{uj)coth[haj/2kBT]cos[ujs]. (6) 



Later we will concentrate on environments that lead to 
periods in time with a negative dephasing rate. 

Recently, we inve stig ated system-environment correla- 
tions of this model |23| and found the useful representa- 
tion 



Ptot(t)= / ^ie-l«l^ 



P{t;tn'»\0{^\ (7) 



of the total state. It represents a partial P-representation 
where the environmental degrees of freedom are ex- 
panded in terms of coherent states |^). Here, ^ = 
(^i;C2,'-') is a vector of complex numbers and we 
consistently make use of the notation d'^^/n := 
d^^i/iT (i^^2/7r- • • (see also jll]). Furthermore, we sym- 
bolically write exp[-|^|2/n]/n := ]\x<^^v[-\^\? Inx]/nx 
involving the mean thermal occupation number fix of the 
A-th environmental mode. The system part of the total 
state is encoded in a matrix-valued partial P-function 
P{t) with values in the 2x2 dimensional state space of 
the qubit. 

In order to represent a solution of the to- 
tal Schrodinger-von-Neumann equation with initial 
Ptot(O) = Psys(O) ® pthorm, the partial P-function in ([7]) 
reads 



-4+(i;C,r)poo 6(t;e,r)poi 1 



Here, = exp[-A{t) ± {{a{t)\^) + i^\a{t))}] and B = 
exp[-int]exp[B{t)-{{b{t)\^) - where we have 

introduced the complex time dependent vectors a{t) = 
{ai{t),a2{t), • • • ) and b{t) with scalar product (a(t)|^) = 
J2x '^xit)^\ and vector components 
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Furthermore, we use the abbreviations 
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Initially, P = psys(O) = /3red(0) and note that there are 
no approximations necessary to achieve the result ([5]) and 
thus via ([7]) to obtain the exact state of the composite 
system (see also [U 



III. RANDOM UNITARY MODELS OF QUBIT 
DECOHERENCE 

With an eye on experimental conditions, decoherence 
of qub its is often modelled by random unitary dynam- 
ics |27l , Isoj . In terms of a dynamical map, this implies 
that there exists a relation 
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with suitably chosen probabilities pk > and unitary 
maps Uk- Indeed, on the level of the reduced state, 
single-qubit decoherence can always be modelled in this 
way 24-26]. It is also worth noting that random uni- 
tary dynamics emerging from an open quantum system 
with environmental initial pure state can always be "un- 
done" (quantum error correction) [10, |4l| . 

The most straightforward random unitary realization 
of single-qubit decoherence with state (U) at time t is 
provided by the simple quantum operation 



1+|PWI ^ ^-^qta. 
l-\Vit)\\ ^_,n,^ 



Prcd(0)e^^*"^ (12) 



which is obviously of the form (jlip employing just two 
unitaries Ui — exp{—iD,taz/'2}, U2 = cxp{—intaz/2}az 
and probabilities pi.2 = (1 ± \T>{t)\)/2. Recall that ac- 
cording to (O, \V{t)\ = exp{— /p 7(s)ds}. 

It is worth noting that the very same formal relation 
holds true for any two-time map £{t,t') : Picd{t') 
Prcd{t) such that 



Prod(^) = 



1 Ui(t-t )pr:cd(t )U{{t - t ) 

U2{t-t')p,,A{t')Ul{t-t') 

(13) 
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with \'D{t,t')\ = exp{— /j*, 7(s)ds}. As 7(5) need not 
be positive for all times (see later) , the prefactor of the 



second contribution. 



2 , may turn negative for t' 
and t near times of negative 7(5). Thus, for such times 
{t\t), the map £{t,t') in the form (fT3|) ceases to take 
the form of a random unitary map. Indeed, using the 
Jamiolkowski isomorphism [i^] it is straightforward to 
see that \V{t,t')\ < 1 or 



7(s)ds > 
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is a necessary and sufficient condition for the dephasing 
map S{t,t') defined above to be completely positive [171] . 

The random unitary form ((T2|) of the dephasing map 
£{t,0) is simple but has the drawback of not represent- 
ing an intuitive dynamical picture of the process due 
to the time dependence of the probabilities pi = pi{t) 
and p2 — P2{i)- Here, we want to develop an alterna- 
tive random unitary representation in a more systematic 
way that opens the door for further generalizations as 
explained below. 

Recall that pure dephasing of a system in the ba- 
sis implies a controlled-unitary form Utot{t) = 
^-iHt^tt/h ^ J2„\n){n\ ® Unit) of the total propa- 
gator ^43. Here, Unit) = e-'""*^'' with i?„ 
(n|i?tot|fi) are system-state-dependent propagators for 
the environment (see also |44|). Pure dephasing 
is then given by the dynamics (?T.|prcd(Ok'^) = 
trcnv(C/„(t)pcnv(0)C/tj(t)) (n | p^od (0) | m) . In our case 
of a single qubit there is a single decoherence factor 
V{t) —tTcnv{Uo{t)pcnv{0)Ul) 'AS in ([5]). The two propa- 
gators are determined by the environment Hamiltonians 
= {i\Ht,t\t) with Hi = hn/2 + J2^^^ hgxial + ax) + 

X^aLi t^xO'\o,x and two sign changes for Hq. We next 
employ the Wigner representation of the environmental 

initial state Wo (a, a*) = / e€*"-«"*tr (e«"^-«*"po„v) 

and for the operator ulit)Uoit), accordingly. The lat- 
ter's Wigner Weyl symbol we denote by J7oi(a, a*) = 



J 7T 



tr 



(^e^-''i'-Ulit)Uoit) 



Wc find 



V{t) = Woia,a*) Uoiia,a*,t). 



(15) 



Due to the harmonic properties of the environment, 
the corresponding propagators Unit) are known ex- 
plicitly and lead to the phase factor Uoiia,a* ,t) — 
exp{— i$(a, a*, t)} with the phase 



'ds + c.c. (16) 



We see that P(t) is just an average over a random com- 
plex number of unit norm. For a thermal initial state the 
initial Wigner function Wq = exp{ — |a|^/(n -I- ^)} is 

positive. Thus, expression (jl5p corresponds to the ran- 
dom unitary representation for Sit, 0): 



Prcdit) = 



Woia, a*) U^it) prcd(O) U^it). (17) 



This corresponds to a random evolution 
Ua = cxp{— « Hais)ds/h} of thc qubit 
with diagonal stochastic Hamiltonian H^it) — 

- Ex 9xiax /o e-'^>-'ds + al /J e^'^^'ds^ . Note 

that in this representation the probability of occurrence 
of a particular unitary evolution is given by the value 
of the initial Wigner distribution and is thus time 
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independent. In this sense, the second random unitary 
representation p7p reflects an ensemble of experiments 
where the system dynamics is influenced by some 
(classical) stochastic process. Note that both random 
unitary representations of the dynamical map £{t, 0) are 
exact - no restriction on the sign of 7(5) is necessary. 

It may appear tempting to define a two-time 
map T{t,t') through ^ with Ua Ua{t,t') = 
exp{— i J^, Ha{s)ds/h}. However, it is clear that 
J'{t,t') ^£it,t') unless t' =0. 

We close this section pointing out an interesting 
additional observation: the random unitary represen- 
tation (jl7l) for the quantum dephasing model is not 
restricted to single-qubit-dephasing. In fact, for an 
arbitrary system Hilbert space dimension, the very 
same construction works for all dephasing factors 
T^nmit) =trenv(t^n(i)/5env(0)J7jti) ^ quautum osciUator 
environment model. So even for larger Hilbert space di- 
mension than two - on a local level - pure dephasing 
based on a quantum oscillator model like ([1]) cannot be 
distinguished from random unitary dynamics. For gen- 
uine quantum decoherence, one needs "more quantum 
mechanical" environments 12611. 



IV. MODEL OF DECOHERENCE WITH 
PERIODS OF NEGATIVE DEPHASING RATE 

The physics of the harmonic oscillator environment 
model is encoded in its spectral density J(i^). As we are 
here interested in instances of negative dephasing rate, 
we chose a particular superomic spectral density with 
sharp cutoff at frequency uJc 



J(lu) — K—Q{ll> — UJc) 



(18) 



and, for simplicity, concentrate on the high-temperature 
limit ksT ^ hwc- In (jl8p . n denotes a dimensionless cou- 
pling constant. In the high-temperature limit the time 
dependent dephasing rate © can easily be evaluated an- 
alytically, we get 

h^[t) = 8.ksT ( ^^"l . (19) 



(uJct) 



As can be seen in Fig. [TJ ^(t) turns negative in certain 
restricted periods, while the integral 



^{s)ds 



sin(a;ct) 
{ujct) 



(20) 



stays positive, as expected for a completely positive map. 

Though many of our results do not rely on any special 
choice of J(w), in the following, whenever we show fig- 
ures, we will use the spectral density ((TH]) and account for 
the high-temperature limit by choosing T = lO/iWc/fcs. 
Furthermore, we choose k = 10~^ troughout this paper. 
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Figure 1: While the dephasing rate 7(f) clearly shows do- 
mains of negative values, the integrated quantity '){s)ds 
stays positive, ensuring complete positivity of the dynamics. 



V. TOTAL STATE AND PURITIES 

Coupling to an environment leads to the built-up of 
correlations between system and environment and thus 
to changes in local entropies. From the point of view of 
information theory, regarding the reduced dynamics as 
a "channel" for quantum information, several quantities 
related to entropy are of interest [43 • Unfortunately, 
since all these quantities (quantum mutual information 
etc.) are based on von Neumann entropy -tr(plogp), 
they are hard to compute, unless one deals with very 
small systems or Gaussian states. 

In a very first step we here want to determine purity 
P=tr(p^) as an indicator for the mixedness of states, 
which is related to the linear entropy via Sl = I — P. 
Clearly, as with entropy, total purity Ptot is preserved 
under unitary evolution with -fftot ■ The sum of local pu- 
rities Psys(i), Penv{t) howcver, will be smaller as t > 0. 
For the initial product state we have Ptot = PsysPenv and 
it appears natural for alH > to consider the difference 
of logarithms of P as a simple measure of "information 
correlation" 



C{t) ^ log(Ptot) - l0g{Psys{t)) - log(Pe„v(t)) 



= log 



= a 



sys 



^sys(O) 
Psys(^) 
(t) + Ce„v(t) 
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^cnv (^) 
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Here the contributions Cgys and C, 



correspond to the 
amount of correlations created between system and envi- 
ronment due to the increase of the local entropies in the 
two subsystems. 

Having the total state ([7]) at hand, all these quantities 
can be determined easily for our dephasing qubit. For 
instance, the qubit purity is readily determined to give 



Psys{t) 



I + + {x^ + y')\V{t)\ 



(22) 
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Figure 2: Dephasing rate ^{t) (solid) and information correla- 
tion C{t) (dashed) against time for a qubit with initial purity 
r — 0.98. While C{t) shows a monotonic increase for posi- 
tive dephasing rates (7(f) > 0), we observe a decrease of the 
information correlation for y{t) < (highlighted domains). 



Here and in the following we denote by r = 
{x,y,z) =ti{ap) the coordinates of the Bloch vector of 
the initial state of the qubit. Somewhat more involved, 
yet still easy to determine is the purity of the environ- 
ment. We find 



PcnAt) = ^ (1 + + (1 - ^^)\G{t)n Pc, 

with the initial environmental purity 



.(0) (23) 



logPcnv(0)= / dujJ{ij)log{taTih{huj/kBT)). 
Jo 



(24) 



In ([22]) , the time dependence arises from the decoherence 
factor — exp(— j{s)ds) of qubit dephasing with 

the rate j{t) — 4 J^ds J^dujJ{uj)coth.[huj/2kBT]cos[uJs] 
from ([5]). By contrast, for the environment 
the time dependence is governed by a fac- 
tor \G{t)\ — exp{— J^T{s)ds) with a dual rate 

r{t) = 4 J^ds J^dujJ{uj)tiinh[fiaj/2kBT] cos[ws]. 

The rate of change of the correlation C{t) from ((2T|) 
stems from the two contributions C — Csvs + Ccnv with 
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where the initial state of the qubit determines the factors 
a = (14- z^)/{l - z2) and 6 = (1 + z^)/ix^ + y^). 

In the high temperature limit, by means of (j25p 
and (j26p all quantities in ([2T|) can be obtained read- 
ily. Reflecting the huge dimension of the environments 



Hilbert space, it turns out that the environmental con- 
tribution is very small compared to Csys(i), and hence, 
we have C{t) « Csys(i). Therefore, from ((25|) one would 
expect C{t) to be correlated with 7(t), independently of 
the special choice of J(a;). In Fig. [5] we compare the 
information correlation to the dephasing rate. Clearly, 
a change in C depends on the sign of 7, and thus C{t) 
decreases in the domains of negative j{t). 



VI. QUANTUM AND CLASSICAL 
SYSTEM-ENVIRONMENT CORRELATIONS 

Having access to the total state we can also investigate 
the nature of system-environment correlations. In earlier 
work we have shown that quantum correlations need not 
exist in such open system models, in particular in the 
high-temperature limit. In such a case, the total state 
may still be separable and thus does not rely on any 
exchange of quantum information. Here we argue very 
much as in |23| . 

With the time and temperature dependent function 



S{T, t)^A f ds f dT f 
Jo Jo Jo 



duj 



Jo Jo (27) 
X J(w) cxp [ftw/fcr] cos[w(s — t)] 

we have shown that the total state is separable, as long 
as 



S{t) < In , 



(28) 



In the high temperature limit, with our special choice of 
J(w) this quantity can be easily evaluated yielding 



S{T,t) =4k 



1 sin[wci] cos[a;c^] 1 

2 UJrt UJrt^ UJrf^ 



(29) 



On the other hand, the total state turns into an en- 
tangled state at the latest when its partial transpose 
yields a negative expectation value (^Ip^^J'lvt^) in some 
state I^E") of the composite system [4a]. By means of the 
representation ([7]) we have shown that with the time and 
temperature dependent function 



EiT, t)= 8 ds dr 



du! 

JO JO (30) 
X J {u) sinh [hoj / kT] cos[uj {s — t)] 



entanglement is proven for 
E{T, t) > In 



y 



(31) 



Therefore E can be interpreted as "entanglement poten- 
tial" , which implies quantum correlations if it exceeds the 
initial state dependent threshold given in (1311) . Again, in 
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Figure 3: E(i) (solid) and Ccnv(i) (dashed) against time for 
a qubit with initial purity r — 0.98. Obviously, a change in 
the local entropy Conv is accompanied by a correlated change 
in the entanglement potential E. Note, that due to the con- 
dition H28|) for separability of the total state, an increase of E 
does not lead to system-environment entanglement. 



First of all, we want to stress that due to the high tem- 
peratures investigated here, the total state of system plus 
environment is separable at all times. This applies even 
if we choose an initial state in the equatorial plane of 
the Bloch sphere with purity of r = 0.98, as can be seen 
from (051) . Nevertheless, it becomes obvious from Fig. [3] 
that the loss of information in the environment, described 
by Ccnv, is correlated with the magnitude of the entan- 
glement potential E. However, these correlations are not 
of quantum nature, as already anticipated. Only if we 
employ even higher purities (r = 0.997) like in Fig. |4j 
condition ([^5]) ceases to hold true and entanglement can 
be detected in certain domains of time. 

Let us stress that the relation between the two consid- 
ered correlation related quantities exists, again, indepen- 
dently of the special choice of J{uj): since in the limit 
ksT ^ hujc the general expressions for the derivative 
E{T,t) of the entanglement potential and the dual rate 
r{T,t) essentially coincide, also the chang es in Ccnv and 
E are correlated. 
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Figure 4: Ccnv(i) against time for a qubit with initial 
purity r = 0.997. The highlighted domains correspond 
to time intervals where the entanglement potential E > 
ln[(r-z2)/(a;2-H2/2)]. Clearly, an increase of the local en- 
tropy Ccnv(i) here leads to the growth of quantum correlations 
between system and environment. 



the high temperature limit the special choice of J(w) en- 
ables us to obtain the following analytical expression 



£;(T,t) = 



I COs[lLJc^] 
3 [i^rtf 



+2 



sin[a;ct] sin[ajci] 



(32) 



VII. CONCLUSIONS 

We have investigated the non-Markovian dynamics of 
a decohering qubit and its environment. Since it can be 
modeled by means of random unitary evolution, we have 
stressed that on the level of reduced dynamics, there is 
no need for a quantum description. 

Nevertheless, considering the full dynamics of sys- 
tem plus environment, we have investigated the measure 
C — Cgys + Ccnv for system-environment correlations that 
emerge from an increase of the local entropies of the two 
subsystems. We have found this quantity to be closely 
related to the dephasing rate ^{t), reflecting the non- 
Markovian character of the dynamics. 

Furthermore, we have found that the environmental 
contribution Conv to the "information correlation" C is 
related to the "entanglement potential" E investigated 
previously in Though in most cases the total state 
of system and environment turns out to be separable, we 
have also shown that, even at high temperatures, quan- 
tum correlations can occur in the total state if the initial 
state of the qubit is sufficiently pure. 

We are confident that our approach will be helpful for 
further investigations with respect to information flow in 
open system dynamics. 



We now use these expressions to relate the build up 
of correlations between system and environment to the 
environmental contribution Ccnv to the information cor- 
relation, introduced in the preceding section. 
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